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ABSTRACT 

The basic and intrinsically disordered C-terminal do- 
main (CTD) of the linker histone (LH) is essential 
for chromatin compaction. However, its conformation 
upon nucleosome binding and its impact on chro- 
matin organization remain unknown. Our mesoscale 
chromatin model with a flexible LH CTD captures a 
dynamic, salt-dependent condensation mechanism 
driven by charge neutralization between the LH and 
linker DNA. Namely, at low salt concentration, CTD 
condenses, but LH only interacts with the nucleo- 
some and one linker DNA, resulting in a semi-open 
nucleosome configuration; at higher salt, LH inter- 
acts with the nucleosome and two linker DNAs, pro- 
moting stem formation and chromatin compaction. 
CTD charge reduction unfolds the domain and de- 
condenses chromatin, a mechanism in consonance 
with reduced counterion screening in vitro and phos- 
phorylated LH in vivo. Divalent ions counteract this 
decondensation effect by maintaining nucleosome 
stems and expelling the CTDs to the fiber exterior. 
Additionally, we explain that the CTD folding depends 
on the chromatin fiber size, and we show that the 
asymmetric structure of the LH globular head is re- 
sponsible for the uneven interaction observed be- 
tween the LH and the linker DNAs. All these mecha- 
nisms may impact epigenetic regulation and higher 
levels of chromatin folding. 

INTRODUCTION 

The genome in eukaryotes is encoded in long double- 
stranded deoxyribonucleic acid (DNA) molecules that are 
densely packed with proteins to form the chromatin fiber. 
Chromatin's hierarchical structures govern transcription, 
replication, and repair of the genome (1,2). At the first level 
of organization, the DNA makes about 1.75 turns around 



the histone octamer (two copies each of the core histones 
H2A, H2B, H3, and H4) to form the nucleosome. Linker 
DNA connects consecutive nucleosomes to form the 'beads- 
on-a-string array', which compacts at physiological salt 
conditions with added linker histones HI /H5. Higher levels 
of organization have remained controversial over more than 
three decades of research (3,4), inviting many experimental 
and computational studies (5). 

Linker histones (LHs) are essential in higher organisms 
(6), regulating the transcription of genes and chromatin ar- 
chitecture (7). Depletion of LHs induces an expansion of 
chromatin in vivo (8,9) and destabilizes chromatin fibers in 
vitro (10,1 1). The compaction effect of LHs is attributed to 
the emergence of nucleosome stem motifs (12,13), which fa- 
vor zigzag folding in 30-nm fibers (1 1,14,15). Because these 
processes depend on the reorganization of the unfolded LH 
regions upon binding to chromatin, an understanding of 
this mechanism is crucial. 

LHs consist of a globular head (GH), and N-terminal 
(NTD) and C-terminal (CTD) domains (6,16). The -80 
amino acid GH was determined by X-ray (17) and nuclear 
magnetic resonance (NMR) (18) for avian and yeast LHs 
to be in the 'winged helix' family of DNA-binding proteins, 
conserved among HI and H5 histones (19). The GH is es- 
sential for the binding of the LH in the nucleosome and sta- 
bilizes 10-20 bp of linker DNA (20). The exact position of 
the GH in the nucleosome is uncertain, although most stud- 
ies indicate that it binds the nucleosome at the dyad axis 
or nearby (i.e. 1-2 DNA turns away) (16,21-23). The 20-35 
amino acid NTD is intrinsically disordered and does not ap- 
pear to affect higher-order chromatin compaction (24). The 
— 100 amino acid CTD is lysine-rich and intrinsically disor- 
dered (25), accounting for most of the sequence heterogene- 
ity among LHs (26). The CTD is crucial for the compaction 
of chromatin (24) and is responsible for the high-affinity 
binding of LH to chromatin in vivo (27). Recent FRET ex- 
periments have shown that DNA binding triggers the con- 
densation of the CTD, with folding dependent on whether 
the LH binds to the nucleosome or to naked DNA (28). In 
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Figure 1. Integrated chromatosome coarse-grained model with refined 
linker histone. (A) The nucleosome core includes the wrapped DNA and 
is modeled by 300 optimized surface pseudo-charges (smallest spheres col- 
ored in pink, magenta, and blue); linker DNA beads (large red spheres, six 
beads for 209bp NRL) are derived from the worm-like chain model; his- 
tone tails H2A1 (yellow), H2A2 (orange), H2B (magenta), H3 (blue), and 
H4 (green) have resolution of five amino acids per bead. The linker his- 
tone GH is represented by six rigid beads (gold) based on the full atomic 
structure of HI. 4 (A and B), and the CTD is modeled as 22 flexible beads 
(cyan) of the same resolution of core histone tails (A and B). (B) Full-atom 
LH model shown as a ribbon representation and by Poisson-Boltzmann 
electrostatic surface (21,29) and our coarse-grained model with flexible 
CTD; the electrostatic representation has red-white-blue scale with range 
[-10e,10e]. (C) Definition of amino acid pairs A and B for the rat HI. 4 
CTD, used in distance measurements in experiments (28). 

addition, distances measured between CTD residues devi- 
ate from those predicted by CTD full-atom models (29). 
These clues suggest that CTD folding is dynamic and de- 
pends on both DNA binding and the chromatin environ- 
ment. 

Computational modeling has become essential for com- 
plementing experimental findings of biomolecular struc- 
tures and associated processes (30,31). Modeling can help 
elucidate the atomic and molecular mechanisms behind the 
structural as well as regulatory roles of LHs in chromatin. 
Various chromatin models have characterized many key el- 
ements that control chromatin structure (reviewed in (5)). 
However, most computational and theoretical models, in- 
cluding our prior model, treat the effects of LH crudely or 
neglect them altogether. In turn, high-resolution LH models 
reveal neither the associated folding dynamics of the CTD 
nor the effect of LH in the global reorganization of chro- 
matin upon LH binding (21,32). 

Here, we enhance our mesoscale chromatin model with 
a refined GH and flexible CTD (see Figure 1). This re- 
fined LH model captures the spontaneous condensation of 
the CTD upon binding the nucleosome and reveals a dy- 
namic charge neutralization mechanism between the CTD 
and linker DNAs. We also observe a synergistic role be- 
tween the LH domains and nucleosome stem formation 
in chromatin fibers and explain the observed uneven or- 
ganization of the linker DNAs by the intrinsic asymme- 
try of the GH. The CTD expansion that we discover for 



smaller oligonucleosomes clarifies size variations observed 
experimentally. In addition, we show that unfolding of 
weakly charged CTDs decondenses chromatin, support- 
ing the regulatory role of LH's post-translational modifi- 
cations in transcription. These insights into the mechanism 
of CTD and chromatin condensation suggest specific ionic 
and amino-acid modifications that control chromatin or- 
ganization and hence genomic accessibility and regulatory 
processes. 

MATERIALS AND METHODS 

Our coarse-grained model has characterized important 
structural mechanisms in chromatin, such as the role of elec- 
trostatics and core histone tails in fiber compaction (33), the 
dependence of chromatin architecture upon nucleosome re- 
peat length (NRL) interfiber (15) and intrafiber variation 
(34), and the compaction and heteromorphicity induced by 
LH and Mg 2+ in equilibrated fibers (11) as well as fibers 
under tension (35,36). Our model also predicted that fibers 
with NRLs close to 177 and 187 bp adopt 30-nm zigzag 
structures (15,37), as recently revealed by high-resolution 
cryo-EM 3D reconstructions (38). 

Our mesoscopic model has been developed for more than 
a decade, and our modeling strategies were recently de- 
tailed in references (15,33). In brief, the nucleosome core 
is approximated as a rigid surface made of 300 pseudo- 
charges, adjusted using Discrete Surface Charge Optimiza- 
tion (DiSCO) at different salt concentrations (39). The core 
histone tails are modeled as flexible charged beads with a 
resolution of five amino acids per bead using the Warshel- 
Levitt united-atom protein model (33). The linker DNA 
is modeled using an elastic worm-like chain of spheri- 
cal beads, where n\> = 6 for NRL = 209 bp. The negative 
charge of DNA is assigned at different salt concentrations 
using the Stigter's procedure (40), as detailed by Schlick 
et al in reference (41). All elements interact with electro- 
static and excluded volume terms as described in reference 
(15). The effects of monovalent ions are captured by adjust- 
ing the Debye-Huckel charges on all chromatin elements at 
different salt concentrations to reproduce atomistic proper- 
ties of the electric field (33). The effect of divalent ions at 
low concentration is approximated by reducing the repul- 
sion and persistence length of the linker DNAs (11). The 
model, however, cannot account for ionic effects between 
LH and nucleosomes at high monovalent salt concentra- 
tions (>200 mM) (42). Accordingly, our study focuses on 
physiological conditions below this ionic strength. 

In our prior LH model, the linker histone HI. 4 (rat) 
was approximated using three rigid charged-beads oriented 
through the nucleosome dyad axis (15). Here, we refine the 
LH model to capture the dynamics and flexibility of the 
C-terminal domain as well as the non-uniform charge dis- 
tribution of the globular head (see Figure 1 and Supple- 
mentary Information for additional details). We derive the 
GH from the full atom structure of the rat HI. 4 (21,29). 
A 6-bead model is obtained by combining the Shape-Based 
Corse Grained (SBCG) method (43,44) with DiSCO (15,39) 
(see Supplementary Tables S.l and S.2). We assume a sym- 
metric positioning of the GH in accord with recent high res- 
olution experiments (12) and orient GH in the dyad by the 
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full atom prediction model for GH and CTD of HI (21,29) 
(see Figure 1 and Supplementary Figure S.l). We coarse- 
grain the highly basic and intrinsically unstructured CTD 
to capture the elasticity and electrostatics of the segment 
at different salt conditions using the same strategy applied 
for the core histone tails (see Supplementary Figure S.l A), 
with five amino acids per bead, so that 22 beads represent 
the 1 1 1 -residue rat HI .4 CTD (109-219). The extended con- 
figuration of the CTD is 25 nm long. To fit the CTD in the 
medium linker DNA fibers, we use a compressed but elon- 
gated CTD starting configuration of length 10 nm consis- 
tent with FRET data (6) (see Supplementary Figure S.1B). 
The section Refined LH model in Supplementary Informa- 
tion contains an extended description of the new coarse- 
grained model, including the interactions and parameters 
of the LH. 

We sample oligonucleosome conformations at constant 
temperature using five Monte Carlo (MC) moves to equili- 
brate and capture the multiscale dynamics of the system: (i) 
pivot, (ii) translation and (iii) rotation of DNA beads and 
nucleosome cores, (iv) regrowth of core histone tails, and (v) 
translation of LH CTD beads (see details in Supplementary 
Information). 

We simulate oligonucleosomes of 3, 6, 12, and 24 residues 
with a uniform NRL of 209bp (six DNA linker beads) at 
temperature of 293 K. This NRL corresponds to differenti- 
ated cells, like chicken erythrocites, where the fibers are fully 
loaded with LHs and the exchange rates are slow (45). Ac- 
cordingly, we model fibers with one LH per nucleosome. 
Different salt concentrations of NaCl (5, 15, 80, and 150 
mM) as well as low concentration of Mg 2+ are examined for 
both our prior (11) and refined LH models. Different LH 
CTD charged states are considered for all salt conditions. 
Systems are equilibrated by 30 to 140 million MC steps 
for different starting configurations (Supplementary Fig- 
ure S.2), which represent the classical zigzag and solenoid 
models for 30-nm chromatin fibers (5,15). We simulate at 
least eight trajectories per starting configuration. The sim- 
ulations are long enough to guarantee the convergence of 
the energy as well as global and local quantities as demon- 
strated previously (15,37). The last 10 million steps are used 
for statistical analysis (with terminal nucleosomes omitted 
for LH properties), storing the fiber properties every 10 000 
steps to guarantee the decorrelation of the averaged struc- 
tures. 

We characterize the structure of chromatin fibers with 
several standard properties described in the Supplementary 
Information and reference (15): sedimentation coefficient 
Cs*20,w) (Eq. S.6), packing ratio (pr) (Eq. S.7), dimer (dd) and 
triplet (td) distances, entry-exit (a) (Supplementary Figure 
S.3) and triple (ta) angles, and core-core interaction pat- 
terns (I(k)) (Eq. S.8). The LH CTD is characterized using 
the end-to-end distance (R e to), radius of gyration (R g ), rela- 
tive shape anisotropy (/c 2 ), and distance between amino acid 
pairs (see Supplementary Information for details). 



A Equilibrium configurations 

5mM 15mM 80 mM 150mM 150mM+Mg J * 




5mM * 15mM 80 mM 150mM 150mM+Mg J * 

Qk OC &* O* 




c#ot 



J|| ^1 "=% i=4| <=^b 



C Internucleosome interaction patterns 

15 mM 




0 2 4 6/c>7 0 2 4 6 *>7 0 2 4 6 k>7 0 2 4 6 k>7 0 2 4 6 k>7 
fr-nucleosome fc-nucleosome /(-nucleosome fr-nucleosome fc-nudeosome 




l:Mffl.'il 



Figure 2. Equilibrium structures and chromatosome conformations for 
12 x 209 bp fibers at increasing salt milieus. (A) Equilibrium configura- 
tion and (B) nucleosome projections, including a characteristic chromato- 
some snapshot as well as accumulated (red) and averaged (blue) projections 
of the linker DNA in the nucleosome top view and the dyad plane (side 
view), along with globular head position (yellow) and average CTD posi- 
tion (cyan). (C) Average core-core frequency interactions (with different 
scales) and (D) a condensation trajectory for CTD end-to-end distances, 
i^ete (see also Supplementary Figure S.14C). 



RESULTS 

LH CTD condensation: CTD neutralization mechanism 

Our simulations from extended CTD conformations of 10 
nm in agreement with FRET data (6) capture a dramatic 
condensation of the CTD within chromatin fibers. For each 
ionic state, Figure 2 shows representative fiber forms, top 
(nucleosome) and side (dyad) views of sampled linker DNA 
beads configurations (plus average linker DNA orienta- 
tions), inter-nucleosome interaction patterns, and LH CTD 
end-to-end distances (i? e te) during fiber equilibration. 
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At low monovalent salt (5 and 15 mM), fibers adopt ex- 
tended zigzag structures (Figure 2A). On average, the linker 
DNAs cross one another in the nucleosome plane as ob- 
served in Cryo-EM experiments (14,46) (Figure 2B). How- 
ever, in the dyad plane, they remain uncrossed and form 
an open angle, which decreases with the salt concentration 
(Figure 2B). The CTD condensation is fast upon nucleo- 
some binding (Figure 2D) and, overall, the chromatosome 
adopts a bidente configuration where the LH interacts with 
the parent nucleosome core and the entry linker DNA (Fig- 
ure 2B and Supplementary Figure S.4). 

At higher monovalent salt (80 and 150 mM), chromatin 
becomes more compact (Figure 2 A and B). This is consis- 
tent with the LH-dependent compaction observed in vitro 
on nucleosome arrays of similar size (10,11). In the com- 
pact fiber, an ordered zigzag pattern emerges with a sharp 
core-core peak interaction for second neighbor cores (i ± 2) 
as well as a secondary peak associated to the interdigitation 
with fifth neighbors (/ ± 5) (Figure 2C). This compaction 
also reduces the linker DNA crossing angle in the nucleo- 
some plane and the linker's range of fluctuations (Figure 
2B). Compared to lower salt, the CTD condensation takes 
longer (more MC steps) (Figure 2D), but now the CTD 
bridges both linker DNAs, establishing a tridente interac- 
tion (nucleosome plus two linkers) and forming stems (Fig- 
ure 2B and Supplementary Figure S.4). 

As Figure 2D shows, the CTD reorganizes at all con- 
centrations, from an extended, R QtQ ~ 10 nm, to a con- 
densed, R QtQ ~ 2-3 nm, conformation. The equilibration 
length depends on the salt concentration, requiring 5, 5, 10, 
and 25 million MC steps to reach condensed conformations 
with R QtQ ~ 2 nm at 5, 15, 80, and 150 mM NaCl, respec- 
tively. Figure 3B and C reveal a low average relative shape 
anisotropy, k 2 ~ 0.1, and a low radius of gyration, Rg ~ 
3 nm, indicative of a compact quasi- spherical shape. The 
independence of the CTD average conformation from the 
monovalent salt concentration suggests a charge neutraliza- 
tion condensation mechanism between the CTD and DNA, 
which we test next. 

Indeed, CTD end-to-end distances for CTDs with neutral 
charge as well as one and two-thirds of the original charge in 
Figure 3A show a CTD decondensation upon charge reduc- 
tion, forming an extended state for a neutral CTD (R QtQ ~ 10 
nm). The neutral and extended CTD adopt a radius of gy- 
ration almost twice the original one and a much more elon- 
gated shape, with k 2 ~ 0.5 (Figure 3B and C). The fully and 
quasi-fully charged CTDs instead adopt similar structural 
properties, and the CTD with one-third of the charge leads 
to an intermediate state. In correspondence with our syner- 
gistic CTD/fiber condensation mechanism, Figures 2 A and 
3F show that the CTD charge reduction also decondenses 
and disrupts the chromatin fibers (see also Supplementary 
Figures S.5, S.6, S.7, and S.8). This unfolding resembles the 
behavior of fibers without LHs (see Supplementary Figure 
S.9). 

Recent FRET experiments have investigated the conden- 
sation of the CTD in nucleosomes for HI (H1.0 from Xeno- 
pus laevis) in vitro by measuring the distance between two 
pairs of residues of the CTD at 50 mM: 4.2 ± 0.2 nm for 
pair A (101-195) and 4.9 ± 0.3 nm for pair B (101-173) 
(28). H1.0 is shorter than HI. 4 by 25 residues (196 versus 
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Figure 3. LH CTD structural properties for 12 x 209 bp fibers for dif- 
ferent LH CTD charged states: neutral (dashed red), 1 /3 charged (dashed 
green), 2/3 charged (dashed blue) and fully charged (solid black). (A) CTD 
condensation dependence on CTD charge illustrated by the average end- 
to-end radius, R ete , overall LHs in a fiber. (B) Radius of gyration, Eq. S.l 1, 
(C) relative shape anisotropy, Eq. S.14, (D) pair A (124-218) amino acid 
distance and (E) pair B (124-198) amino acid distance. The error bars cor- 
respond to the standard deviation. In (D and E) we also plot the amino 
acid pair distances measured experimentally with FRET for H1(0) at 50 
mM NaCl (28): pair A exp = 4.2 ± 0.2 nm (D) and pair B exp = 4.9±0.3 nm 
(E). At the bottom of the figure we show equilibrium snapshots at different 
salt conditions for uncharged LH CTD (F). 



219 amino acids). We associate pair A in our model with 
region 124-218 (95 residues starting from the penultimate 
amino acid) and pair B with region 124-198 (75 residues 
from the same starting point of pair A). Excellent agree- 
ment with experiments can be noted for Figure 3D and 
E, which plot equilibrium distances between corresponding 
pair beads, extrapolating values for 50 mM. Interestingly, 
our model and the experiments agree that, for fully charged 
CTDs, pair B is ~1 nm larger than pair A, despite being 
22 residues shorter than pair A (see Figures 3D, B and 4B). 
H1.0 and HI. 4 are both rich in lysines and have an identity 
of 46.4% for pair A (47). We also find that the CTD charge 
reduction increases distances for both pairs A and B, con- 
sistent with the unfolding of the CTD at these conditions. 
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Figure 4. LH CTD structural properties for fibers with 3-, 6-, 12-, and 24- 
oligonucleosomes at 150 mM NaCl. (A) Radius of gyration, Eq. S.l 1, and 
relative shape anisotropy, Eq. S.14, as well as (B) distances between amino 
acid pairs A, G124-K218, and B, G124-A198 (as defined in Figure 1C). 
(C) Typical fiber configurations. 



LH CTD configuration dependence on oligonucleosome size 

Several studies have recently presented varying outcomes 
concerning the organization of nucleosome stems and LH 
CTD conformations using Cryo-EM, hydroxyl radical foot- 
printing (12) and FRET (28) for mono-, di- and trinucleo- 
somes. 

To help interpret observations, we examine fiber proper- 
ties for 3, 6, 12, and 24 nucleosomes at 150 mM NaCl in 
Figure 4, Supplementary Figures S.10 and S.ll. For small 
fibers, we see that the average dimer distance is shorter 
and the triplet distance, triplet angle, as well as entry- 
exit angle are all larger, although the values stabilize for 
12-oligonucleosomes. Figure 4A shows that the structural 
properties of the LH CTD exhibit a more complex behav- 
ior. While the shape anisotropy is small for all oligonucleo- 
somes, average values and associated fluctuations increase 
with fiber size, indicating a richer combination of quasi- 
spherical and slightly elongated CTD conformations. The 
radii of gyration strongly increase for small oligonucleo- 
somes and, for trimers, adopt a value of ~5 nm with a quasi- 
spherical CTD. Despite the divergence in R g , Figure 4B 
shows that the distances between the pairs separated by 95 
(pair A) and 75 (pair B) residues remain relatively constant, 
as found experimentally. More detailed insights into these 
structural changes in CTD folding demand both higher res- 
olution FRET experiments and all-atom modeling. 

That both local structural properties of chromatin fibers 
and the conformation of the LH CTD depend on the array 
size is consistent with our charge neutralization mechanism. 
Our quasi- spherical CTDs for small oligonucleosomes and 
more elongated configurations for larger fibers can be ex- 
plained by greater exposure to solvent for small oligonucle- 
osomes. At the same time, the physico-chemical properties 
of the CTD stabilize interactions between specific residues 



independent of size. We also find that 12-oligonucleosome 
fibers are representative for evaluating structural properties. 

Chromatosome asymmetry and heteromorphicity 

The LH globular head provides additional protection 
against enzymatic and chemical hydrolysis of both linker 
DNAs at the entry /exit region of the nucleosome (12,20). 
Earlier structural experiments suggested that this protec- 
tion is symmetrical, about 10 bp per linker DNA (12), but 
recent biochemical experiments (48) using a refined nucle- 
ase technique account for DNA sequence bias and demon- 
strate that HI confers an uneven protection of 11-12 and 
8-10 bp to the entry and exit linker DNAs, respectively. 
Complementary high-resolution models and experiments, 
however, suggest that in this case the GH occupies a sym- 
metric position in the nucleosome, contacting the dyad axis 
(12,20,22). Is there a source of asymmetry in the GH's pro- 
tection of linker DNAs and, if so, what is its origin? 

Our chromatosome configurations for 12-unit oligonu- 
cleosomes at different salt concentrations in Figure 2B show 
that the average position of the linker DNAs is asymmet- 
ric with respect to the dyad axis near the GH for all condi- 
tions. In the nucleosome (top view) plane, the linker DNAs 
cross at a point that deviates from the dyad axis; in the dyad 
plane (side view), the linker DNAs do not cross and both 
tend to be oriented toward the top side of the nucleosome. 
Thus, the uneven protection provided by the LH is captured 
by our model. To determine the mechanism responsible for 
this asymmetry, at 150 mM NaCl we compare our refined 
model for both the GH and CTD to other LH models: our 
prior symmetrical LH model made of three rigid beads (one 
for GH and two for CTD) oriented through the dyad axis 
(Supplementary Figure S. 1) and a hybrid model of one bead 
for GH with our refined CTD (22 beads). We see from the 
chromatosome projections in Supplementary Figure S.12 
that the two simpler LH models (with a one-bead GH) lead 
to nucleosomes where the linker DNAs adopt symmetrical 
configurations on average in both the nucleosome and dyad 
planes, in contrast to our refined GH model (see Figure 
2B). These findings pinpoint the GH asymmetrical struc- 
ture and charge distribution (see Figure IB) as responsible 
for the asymmetric organization of the chromatosome. Ac- 
cordingly, the asymmetric GH is likely responsible for the 
uneven protection of the linker DNAs against enzymatic di- 
gestion. 

The above biochemical studies are bulk experiments that 
provide an average picture of the chromatosome, but do not 
reveal actual configurations of individual nucleosomes. Our 
average chromatosome geometries are consistent with ex- 
perimental data and offer additional insights into the struc- 
ture of individual nucleosomes in a fiber. The red regions 
in Figure 2B show that chromatosomes span wide range 
of configurations. In fact, even at 1 50mM NaCl, when the 
zigzag pattern is prominent (i ± 2 peaks in the internucleo- 
some patterns), the accommodation of chromatosomes in a 
fiber is non-uniform. Representative individual chromato- 
somes in Supplementary Figure S.l 3 show qualitative devi- 
ations from average states. This non-uniform organization 
of chromatosomes is a general feature of chromatin organi- 
zation that emerges independently on the type of LH model 
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Figure 5. Structural properties for 1 2 x 209 bp as a function of salt for re- 
fined (solid black) and prior (dashed brown) LH. (A) Sedimentation coeffi- 
cient, (B) packing ratio, (C) entry-exit angle (see definition in S.3), and (D) 
triplet angle. Error bars correspond to standard deviation. Experimental 
values for the sedimentation (A) are obtained from 12-oligonucleosomes 
made of 207 bp (1 1) (red) and 208 bp (49) (green), native chromatin from 
rat liver nuclei with an NRL ~ 200 bp (50) (blue), and 12 x 207 bp fibers 
(gray). Other experimental dots are associated with native chromatin fibers 
of chicken erythrocytes and are extracted from references (14) (B and C) 
and (46) (D). (E) Compare internucleosome interactions for 150 mM (top) 
and 1 50 mM with Mg 2+ (bottom) from zigzag and solenoid starting con- 
figurations (six plots) and measured patterns from EMANIC (11). 



used. Such an intrinsic heteromorphic mechanism (5,11) is 
likely enhanced by sequence specificity and regulatory fac- 
tors. 



Chromatin compaction dependence on salt concentration 

Many experiments in vitro utilize fibers of 12 nucleosomes to 
elucidate the structural properties of chromatin (11,49,50). 
Accordingly, in Figure 5A we compare sedimentation co- 
efficients obtained at different salt concentrations to those 
obtained with our mesoscale model. The NRLs for the ex- 
perimental fibers range from 200 to 209 bp, similar to the 
NRL of 30-nm model fibers of chicken erythrocyte in dif- 
ferentiated cells (51). All fibers adopt similar sedimentation 
coefficients that increase with the salt concentration and 
are in good quantitative agreement with our 12 x 209 bp 
mesoscale chromatin model. 

High-resolution structural techniques, like Cryo-EM, 
also characterize the packing ratio of chromatin fibers and 
show the increased density of nucleosomes (through the 
longitudinal axis of the fiber). Our comparison in Figure 5B 
shows excellent quantitative agreement. Figure 2A displays 
the corresponding equilibrium configurations of fibers at 
the different salt concentrations. 

An experimental characterization of fiber local struc- 
tural properties is challenging due to the high compaction 
adopted by chromatin at physiological conditions, but mea- 



surements can be reconstructed from Cryo-EM data (14). 
Figure 5C compares measured experimentally entry-exit 
angles, which decrease upon salt concentration, with our 
results for 12 x 209bp models. Excellent agreement can be 
noted. The reduction of this angle corresponds to the for- 
mation of nucleosome stems. This stem is favored by the 
contact of the GH and CTD with both linker DNAs. (Note 
that our computed triplet angle is very similar to the entry- 
exit angle (see Figure 5D); this useful relationship is signif- 
icant because the definition of the triplet angle is more pre- 
cise than the definition of the entry-exit angle, see Supple- 
mentary Figure S.3). Figures 5C and D show that our com- 
puted values are also in good agreement with recent mea- 
surements of the triplet angle of chicken erythrocyte fibers 
at 5 mM of monovalent salt by Scheffer et al (46) and entry- 
exit angle derived from Bednar et al (14). Dimer distances 
measured at 5 mM monovalent salt are also in good agree- 
ment (see Supplementary Figure S.14A). 

Figure 5 also displays structural properties for the prior 
LH model, which is generally in good agreement with exper- 
iments. However, the refined LH model reproduces better 
most structural experimental data, including the increase of 
i ± 1 interactions when starting from solenoid configura- 
tions, in accordance with EMANIC experiments (11) (Fig- 
ure 5E). 



Effect ofMg 2+ 

Using our simplified treatment of moderate Mg 2+ concen- 
tration (reduction of the DNA persistence length and en- 
hanced interactions between linker DNAs) (11), Figure 5 A 
and B show that divalent ions lead to higher sedimentation 
and packing ratio values, even for uncharged CTDs (see 
Figure 3F and Supplementary Figure S.7). The triplet angle 
and the pattern of the zigzag fiber are similar to the case of 
monovalent salt at 150 mM, although long-range interac- 
tions are stronger. 

Importantly, divalent ions affect the organization of the 
CTD in several ways. For highly charged CTDs, the ra- 
dius of gyration remains stable, but the anisotropy increases, 
leading to a more elongated conformation (see Figure 3C). 
This reconfigured CTD structure has an opposite effect to 
pairs A and B, so that distances become similar, d V A — 
d P B- Our snapshots in Figure 3F and Supplementary Fig- 
ure S.5 show that, at reduced CTD charges, the LH CTD 
escapes from the interior of the highly condensed fibers, be- 
cause the charge neutralization mechanism cannot counter- 
act the steric effects, and entropy favors these highly dy- 
namic states. This important effect points to a potential 
regulatory mechanism in chromatin and shows that Mg 2+ - 
induced linker DNA charge neutralization can mimic the 
stabilization role of the CTD on nucleosome stems. 



DISCUSSION 

Our enhanced chromatin mesoscopic model with flexible 
LHs reveals, for the first time, the synergistic condensa- 
tion between LH CTD and the overall chromatin fiber or- 
ganization. The CTD condensation in the nucleosome oc- 
curs dynamically with other nucleosome components and 
is controlled, in part, by a charge neutralization mechanism 
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with the linker DNAs driven by 30-50 positively charged 
CTD residues. This mechanism explains similar general be- 
havior for different linker histone variants with similarly 
charged CTDs (Hlc, Hie, HI 0 ) (9). Highly charged CTDs 
condense and compact chromatin more tightly, therefore 
favoring genome silencing, consistent with structural data 
and regulation activity reported for the highly charged H5 
in chicken erythrocyte (52,53). Analogously, lower-charged 
CTDs adopt variable positions that depart from the inte- 
rior of compact fibers, enhancing genomic accessibility and 
generating potential long-distance interactions with other 
protein remodeling or nuclear elements. 

A dynamic CTD works in tandem with the GH to sta- 
bilize nucleosome stems by bridging the two DNA linkers, 
leading to a tridente interaction of the LH in the nucleo- 
some (two DNA linkers plus the nucleosome core). Thus, it 
is not a stepwise process whereby CTD folds and then com- 
pacts chromatin structure but rather a dynamic interplay 
between CTD condensation and chromatin compaction. 
For small oligonucleosomes, for instance, the CTD expands, 
which could explain observed extended configurations of 
the CTD upon binding protein-free DNA (28). CTD con- 
formations also depend on its charge, which might promote 
variations in folding observed for different LH types (6,29). 
Interestingly, despite our LH positioning in the dyad axis, 
the tridente interaction with the nucleosome leads to an 
asymmetric chromatosome organization, where the entry 
linker DNA is closer to the LH than the exit linker DNA. 
This is due to the asymmetric structure of the GH and could 
explain recent experiments that measured a related uneven 
protection of the linker DNAs by LH (48). The fact that 
GH favors the interaction with one linker DNA is also con- 
sistent with in vivo mapping of the HI interaction surface 
in the nucleosome (54) and NMR in vitro experiments of 
the chromatosome (23); the latter work also observes that 
one of the two H2A core histone tails is involved in the 
LH-nucleosome complex, in agreement with our model (see 
Supplementary Figure S.15). In addition, Supplementary 
Figure S.15 shows that, besides H2A2 tails, H3 tails ex- 
hibit a significant interaction with parent DNA linkers in 
the asymmetric tridente configuration. 

The net positive charge of the CTD varies among LH 
variants and can be decreased by post-translational mod- 
ifications such as lysine acetylation (55), serine, threonine 
and thyrosine phosphorylation (56), and arginine citrulli- 
nation (57). Our findings suggest that such charge reduc- 
tion, along with acetylation of the lysine-rich core histone 
tails, should enhance more open structures and nucleosome 
accessibility, thus favoring transcription. This charge de- 
pendent mechanism could also explain experiments of LH 
phosphorylation, which mimic the effect of LH unbinding 
in chromatin decondensation (56); it also supports the con- 
densation of chromatin observed in the dephosphorylation 
of HI in differentiated cells (58). Additionally, the impact 
of the CTD net charge reduction is compatible with the de- 
condensation of chromatin observed in vitro in the absence 
of LHs (3). The integrated view of the LH dynamic struc- 
ture and chromatin architecture presented here establishes 
a framework for further interrogation of the structural role 
and regulatory mechanisms of LH variants as well as LH 
modifications in chromatin structure and function. 



SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
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